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Abstract 

The exact value of the Lyapunov exponents for the random matrix product P/v = 
AnAn_i ■ ■ ■ A\ with each A{ = E 1 / 2 ^, where £ is a fixed d x d positive definite 
matrix and a d x d complex Gaussian matrix with entries standard complex 
normals, are calculated. Also obtained is an exact expression for the sum of the 
Lyapunov exponents in both the complex and real cases, and the Lyapunov exponents 
for diffusing complex matrices. 

1 Introduction 

Presently there is a great deal of interest, from the viewpoints of probability theory and 
applications to physics and communications engineering, in the statistical properties of 
large random matrices (see e.g. the recent texts [D [221 HZl 129]). Looking back to the mid 
1980's it would have been fair to say that the same applied to the topic of products of d x d 
random matrices. Thus it was in that era that the foundational probabilistic works of 
Kesten, Furstenberg, Oseledec and others in the 1960's and 70's had matured to the extent 
that a book on the subject was written [I]; that a summer research conference was held with 
this topic dominating the subsequent proceedings [7]; and that a number of applications 
to physics were in the full swing of investigation, culminating in the appearance of the 
research monograph [TUj . 

There have been some present day works that have aimed to combine the contemporary 
interest in the eigenvalues of large random matrices with the topic of products of random 
matrices, by studying eigenvalue distributions of products of random matrices, in the limit 
that the size of the matrices is large El [23j [HI [2Ql [28] - There has also been a good 
deal of present day activity relating to the numerical computation of Lyapunov exponents 
ED El [ED El EDI 6]. Regarding the latter, let 

P N = A N A N - 1 ---A 1 , (1.1) 
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where each A, is a d x d independent, identically distributed random matrix such that 
the diagonal elements of A^A have finite second moments. According to the multiplicative 
ergodic theorem of Oseledec [2TJ [23] , one has that the limiting matrix 

V d := lim (PlP N ) 1 '^ (1.2) 

N— >oo 

is well defined, with d positive real eigenvalues e Ml > e M2 > ••• > e Md . The {//j} are 
referred to as the Lyapunov exponents. Already in 1973 Kingman had nominated methods 
to compute the largest Lyapunov exponent [TS] as an outstanding problem in the field. 
The work [21] solves this problem in the case when each Aj is chosen from a finite set, and 
has positive entries. 

Implicit in the need for efficient computational methods is that it is generally not 
possible to compute the Lyapunov exponents analytically. Some noteworthy exceptions 
occur in the case d = 2; see e.g. [161 [171 19] and references therein. For general d, apart from 
the case of diagonal matrices, it seems that the only exact computation of the Lyapunov 
exponents recorded in the literature is when the are real Gaussian matrices with entries 
independent standard real normals. Then it is a result of Newman [19] that 

^ = I(tog2 + *(£^±i)) (i = l,...,d), (1.3) 

where ^>(x) denotes the digamma function. 

It is the aim of this paper to extend available exact results on the evaluation of Lyapunov 
exponents for certain d x d random matrices. In particular, a closed form evaluation of 
{/ij} is obtained for 

At = £ 1/2 G?, (1.4) 

where E is a fixed d x d positive definite matrix and G1 a d x d complex Gaussian matrix 
with entries standard complex normals. Partial results are also obtained relating to the 
exact evaluation of {/Xj} for Aj = E 1 / 2 ^, where G\ denotes a d x d real Gaussian with 
independent standard normals as entries. The latter supplements the E = 1^ result (11. 3ft . 
and the d = 2 result 

= ~\l + \ log (^Tr E + v/detE) , (1.5) 

where 7 denotes Euler's constant. 

Our first result is the analogue of (ll.3p for complex Gaussian matrices with entries 
independent standard complex normals. 

Proposition 1. Consider the matrix 'product ( ti.il) . with each Ai given by ( |i.^[ ) with E = 7^. 
We have 

p. = -qr(d-i + 1). (1.6) 
This result can be generalized by allowing for general E in (11.41) . 
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Proposition 2. Consider the matrix product U.l\) . with each given by ljl-4\ ) for general 
positive definite S. Let the eigenvalues of be denoted {yj}j=i,...,d- We have 



2 Yli<i<j<d(yj 



det 



J J 



i=l,...,k— 1 
=1, 



[(log % )^ r ] i=1 ,... )(i 



i=fc+l,...,d 
3=l,...,d 



Corollary 1. In the setting of Proposition® one has the sum rule 



1 d 

fJ-i + A*2 H h ^ = - ^ ( - log y m + *(m) ) . 



m=l 



We remark that for d large, the fact that ^{d) ~ logrf tells us that the eigenvalues y m 
must have the scaled form y m /d = Y(m/d), with Y(0) = a PDF on (0, 1), for (11.81) to 
have a well defined average as d — > oo. 

The above results will be proved in the next section. It will furthermore be showed that 
Corollary [T] can be proved independent of Proposition [21 and this will allow an analogue of 
(11. 8p in the case of real random matrices A4 = E 1//2 GJ to be derived (see eq. (12.251) below). 
We give too, in eq. (I2.28p . the evaluation of the so-called generalized maximum Lyapunov 
exponent [10] 

1 



L{q) 



lim — lo£ 



P 



N I 



in the case of (jl.4p . In section 3 we calculate the Lyapunov exponents for 



.4; 



lim 

m — Vrvi J- J- 



cfV™ 1/2 



e •? 



d) 



H x + iH 2 , 



1.9) 



(1.10) 



where ifi and H 2 are Hermitian matrices distributed with density function proportional 
to exp(— TrHf/2wi) and exp(— Trif|/2u> 2 ) respectively. Section 4 discusses features of jii 
as given by ( II. 7p . 



2 Proofs 

2.1 Background theory 

A fundamental characterization of the Lyapunov exponents defined below (11. 2p is that they 
satisfy [2H [25] 

//i + --- + /x fc = sup lim —\ogVo\ k { yi (N),...,y k (N)} (k = l,...,d). (2.1) 

In ( 12. ip yj{N) := P/v?/j(0), the supremum is over all sets of linearly independent vectors 
{l/i(0), . . . , 2/fe(0)} and Vol*,- refers to the (generalized) volume of the parallelogram gener- 
ated by the given set of k vectors. In regards to the latter, with 

B N := [ yi (N)y 2 (N) ■■■y k (N)}, 
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so that Bn is the d x k matrix with its columns given by the k vectors {yj(N)}, we have 

Vol k { yi (N), . . .,y k (N)} = det^jv^) 1 / 2 = det{BlP ] N P N B,) 1 ' 2 (2.2) 

Following [HI [19] the basic fact that makes the computation of (12. ip tractable for ma- 
trices ( II ,4p is that the distribution of the random vector G^x/\x\, x ^ 0, is independent 
of x. Thus with (yi(0), . . . ,y k (0)} a set of k linearly independent unit vectors,, and E^xk 
denoting the dxk matrix with l's in the diagonal positions of the k rows, and 0's elsewhere, 
we have 

N N 

BlP^PsB = J] E T dxk Gf YGZEm = II G %^ G l^ ( 2 - 3 ) 

3=1 3=1 

where G c - k denotes the dxk matrix formed by the first k columns of G c j. Substituting 
in (12.21) . then substituting the result in (12. ip . we see firstly that there is no longer any 
dependence on {?/i(0), . . . , y k (0)}, so the sup operation in (12.11) is redundant. We are then 
left with the expression 

* + ■ ■ ■ + ^ = & e n log det ( G ii EG ^) • (2 - 4) 

3=1 

But each Gj independently belongs to the set of complex rectangular dxk Gaussian 
matrices A/J xfc (0, 1) in which each entry is a standard complex normal. The law of large 
numbers tells us that the limit in (12 ,4p can be evaluated as an average over this set, 

/ i 1 + --- + /i fc = (logdet(G^SG c fc ) 1/2 )^ (2.5) 
2.2 Proof of Proposition [I] and Corollary [I] 

Proposition [T] relates to the case £ = I d . Now the set of matrices A/^ xfc (0, 1) have proba- 
bility density function proportional to e~ Tr( - G j tG A Thus the average in (12 .5^ is a function 
of G^G c j. Introducing the complex Wishart matrix W = G c ^G c p we know that the corre- 
sponding Jacobian is proportional to (det A) d ~ k (see e.g. [121 Eq. (3.23)]). Making use too 
of the simple identity 

■I- (det Wy i = Tr log det W (2.6) 
d[i fi=o 

it is therefore possible to rewrite ( 12. 4p in the case £ = Id as 

Hi + • • • + to = ^(( det WT +d ~ k e-™)^ (2.7) 



w>o 



where the average is over all positive definite k x k complex Hermitian matrices. 
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We see that (12. 7p is a function only of the eigenvalues of W. Changing variables to the 
eigenvalues and eigenvectors (see e.g. [HI Proposition 1.3.4]) gives 

Id 

/ii H V^k — tt^ -i-Z d>d - k+IJi , (2.8) 

where 

/oo /*oo ^ 

dxi • • ■ / tted J] II " x ^ 2 ' ( 2 ' 9 ) 

But the integral Z C) & is a particular limiting case of the Selberg integral and as such has 
a product of gamma function evaluation (see e.g. [121 Prop. 4.7.3 with f3 = 2]), telling us 
that 

i d k ^r(d-k + v + i + j) 

^ + - + ^= 0.77.11 



2dfxf± o T(d-k + l+j) 

1 fc-i 1 fe-i 

= -£>(d-fc + l+j) = -£>(d-j). (2.10) 

3=0 j=0 

The result ( II. 6p is now immediate. 

Let's now consider Corollary [TJ We thus want to evaluate ( 12. 5 p in the case k = d 
but with X a general positive definite matrix. For d = k we can write det(G^EG^) 1 / 2 = 
det(G£G^E) 1//2 since then G c k G c k has full rank. Furthermore the probability density corre- 
sponding to A/J xfc (0, 1) is then proportional to exp(— Tr G c k G c k ). Introducing the complex 
Wishart matrix W = G c d G c d \ for which the corresponding Jacobian is just a constant, we 
therefore have 

»! + ■■■ + Hd = ^(logdet£ + (log(detWOe- w \ Y 

2 V \ / W>0/ 

Noting that the latter average is just (I2.10p in the case k = d gives (jl.8p . 

2.3 Proof of Proposition [2] 

Writing X = E 1 / 2 ^, (ESJ) can be rewritten 

/ii + --- + /i fc = (logdet(xtx) 1/2 \ . (2.11) 

Now the set of matrices S _1 / 2 A/^ xfc (0, 1) have probability density function proportional 
to exp(— Tr(XX^S -1 )). The d x d matrix XX' has rank k < d, and so has d — k 
zero eigenvalues. This feature, in the case d < k, distinguishes the set of matrices 
S _1 / 2 A/^ xfc (0, 1) from A/J xfe (0, l)!!^ 2 . Thus the latter has probability density propor- 
tional to exp(— Tr(X^ XY,^ 1 )) with X^X of full rank and thus having no zero eigenvalue. 
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These matrices are termed complex Wishart matrices of mean zero and covariance 
while the matrices XX^ are sometimes referred to as complex pseudo Wishart matrices of 
mean zero and covariance E. 

Let p(i)(A;S) denote the eigenvalue density of the nonzero eigenvalues {\j}j=i,...,k f° r 
the ensemble of complex pseudo Wishart matrices so specified. Noting that 

logdet(xtx) =-^logA, (2.12) 

i=i 

and so is a linear statistic in {Aj}, it follows from ( 12. lip 

1 f°° 

IH + --- + li k =- (logA)p ( i)(A;E)dA. (2.13) 
^ Jo 

Studies in wireless communications (261 123 H3] have required the same averaged linear 
statistic, generalized so that log A log(A — z). In these references, methods involving 
integration over the unitary group have been used to find an explicit expression for this 
generalization of the RHS of (I2.13P in terms of a d x d matrix. Thus with the eigenvalues 
of S -1 denoted as in Proposition [2] it is shown 

poo 

/ log(A - *)p (1) (A; S) dX = C M (M) det t L ™] ( 2 - 14 ) 

J m=l 

where L m is the d x d matrix with entries 



/ °° log(t - z)e~^H k - i dt, i 



m 



and 



{L m ) ij ={ yp«-^>r(k-i + l), i^m,i<k (2.15) 
yj-\ i > k 

(_l)k((k+l)/2-d) Tjd k 

CkAivj}) = 7^r-= L f^- V (2.16) 

1L=i tl ([h< t<j <d(yi -Vj)) 

For purposes of computing (11.91) . and also for purpose of making the presentation more 
self contained, it is of interest to revise the derivation of (I2.14p . The first step is to consider 
d x k matrices X with a probability density function 

Changing variables to the nonzero eigenvalues and eigenvectors of XX' by writing XX' = 
?7diag (Ai, . . . Afc, 0, . . . , 0)17^ for U unitary shows 

Pfe(Al ' • • • ' AfcH nSr (2+ w-Hi + i) n ^ nj^ - a,) 2 

x lim Afc+liA ^ J exp(-Tr(t/diag(A 1 , . . . A^E" 1 )) (WdU), (2.17) 
6 



where (U^dU) denotes the Haar volume form for dx d unitary matrices normalized so that 
J (WdU) = 1. The matrix integral is precisely what is known as the Harish-Chandra- 
Itzykson-Zuber integral (see e.g. [121 Prop. 11.6.1]), which can be evaluated as a determi- 
nant to give 



det£- fc n?=or(l + /) 



Pfc(Ai, . . . , Afc) 



nto 1 r(2 + i)r(d -k + i + i)t\ 1£f<lst 



(vj - yi) 

A fe+ i,A d -)-0 ^ j^-j;^ , (Aj — A;) 



x lim 

l<j<K<f 



detfe* v 



i,i=i,...,d- 



(2.18) 



The limit can be carried out by power series expanding columns k + 1, . . . , d, leaving us 
with the explicit determinant formula 



P fc (Ai, . . . , Afc) 



■ n (^-v 



x det 



i<j<i<k 



i=l,...,k 
j=l,...,d 



i=d+l,...,fc 

i=l,...,d 



(2.19) 



The second step is to use (I2.19P to compute the average 



(2.20) 



For this, one notes that Pk consists of two anti-symmetric factors in {Xj}. Since, according 
to the Vandermonde determinant formula 



n (A i -A I )=^n A ^ 

l<j<l<k 1=1 

where A denotes the anti-symmetrization operation, we can replace Yli<j<i<k(^j ~ ^i) * n 
the integrand implied by (12.201) by k\ nf=i The integrals over {\j} can now be done 
row-by- row in the remaining determinant, and we obtain 



IK A ; 



p k nt/ i] rTi<j<i<d(2/i - yi) 



x det 



[Vj l ]i=k+l,-,d 
J j=l,...,d 



(2.21) 



The third and final step is to differentiate this formula with respect to fi and set 
/i = 0. On the LHS this gives the LHS of (I2.14p . Recalling that the differentiation of a 
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determinant with respect to a parameter is equal to the sum of determinants with a single 
in each differentiated, we see that the RHS of ( 12.141) indeed follows by performing this 
operation on the RHS of (12.211) . 

Setting z = in ( 12. 14|) allows the integral in row i = m to be evaluated. Doing this 
and also taking out appropriate common factors from each of the first k rows shows 



logAp(i)(A; S) dX 
LiU((W)M 

E det 



T\-l<i<j<d(yi Vj) m= i 



i-ll 



i=l. . . . , m — 1 



[2/ 



d-\-k—i 



i — m+1 , . . . ,d 
*=l,...,d 



.(2.22) 



Furthermore, reversing the order of the rows i — k + 1, . . . , d and recalling that in general 
a determinant with a single row having each entry a sum of two terms is equal to the sum 
of two determinants, it follows from (I2.23P that 



logA P(1) (A;S)rfA = rf - 

i-ll 



E 



- det 



i — 1 , . . . ,m — 1 



i = 7Tl + l , . . . ,( 

3=1,. ..,d 



+ V(k - m + 1) det^f 1 ]^^!,...^ . 



(2.23) 



Substituting in ( 12. 13f) and making use of the Vandermonde determinant evaluation gives a 
result equivalent to ( 11. 7ft . 



2.4 Second proof of Corollary [T] 

Comparison of (ll.7p and (II. 8ft shows that it suffices to check that 



n 



l<i<j<< 



i(Vi - Vj. 



E det 



m=l 



i-ll 



[y 



i-ll 



i — 1 , . . . ,m— 1 
3=1 d 
m—ll 

j J I 

I i — m+ 1 1 ■ ■ ■ ,d 



3=1, ...,d 



m=l 



(2.24) 



For this purpose, we observe that the sum over determinants on the LHS is equal to 

d 



But this expression in turn can be rewritten 

d 



fi=0 



.7=1 



j=l m=l 

The identity (I2.24p now follows from the Vandermonde determinant evaluation. 
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2.5 The sum rule (ECS) in the case of A { = Y}I 2 G\ 

Suppose the A» in (II. ip are given by Ai = Y}I 2 G\, where each G\ denotes a d x d real 
Gaussian with independent standard normals as entries. Let A/"J xd (0, 1) denote the set 
of real rectangular d x k Gaussian matrices in which each entry is a standard normal. 
Analogous to (12. 5 p we have 



+ ■•• + /** = (log det (G r k T XGl 



1/2 



^ r xfe (o,i) 



where G\ denotes G\ restricted to the first k columns. As in the argument below ( 12. 10[) in 
the complex case, for d = k we can write det (G k T T,G r k Y^ 2 = det(G T k G r k 'Z) 1 ^ 2 . Furthermore 
the probability density corresponding to A/"J xfc (0, 1) is then proportional to exp(— Tr G T k G r k T /2). 
Introducing the real Wishart matrix W = G T d G d T , for which the corresponding Jacobian 
is proportional to (det W)' 1 ^ 2 , we therefore have 



fii + • • • + = 2 (log det E + ((det W)~ 1/2 log(det W) e" Tr 



W/2 



W>0 



where the average is over positive definite k x k real symmetric matrices. To evaluate this 
average we make use of (12. 6 p to write 



log(det W) e- TrW ' 2 ) = ^-/(det WY'^e^/ 2 

I w>o dfi\ i w>o 

Arguing now as in the derivation of ( 12.81) shows that the RHS is equal to 

1 d a 



Zd -i/2 



J d -1/2+/1 



H=0 



where 



dxi 



n 

j=l l<j<l<d 

Like (I2.9p . this integral has an evaluation in terms of a product of gamma functions 
(see e.g. p21 Prop. 4.7.3 with (3 = 1]), giving 



log(det W) e 



-TrW/2 



w>o d[x 



d_ d /jl I> + l/2 + i/2) 

^ 11 r(i/2 + J /2) 



M=0 



d-1 



= dlo g 2 + £>((; + l)/2). 
Consequently, the sought analogue of (I1.8P is 

j d-1 

\ii + ■ ■ ■ + lid = -= ( log det E + d log 2 + *(C? + i)/ 2 ) 



(2.25) 



j"=o 
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2.6 The generalized maximum Lyapunov exponent 

For the matrix norm in (11. 9ft we take ||Pjv|| = su P^ : |^|=i \Pnx\- This gives 

L(q) = sup S:]Sl=1 log (\P N x\ q y 
Arguing as in the derivation of (12.51) . for the A{ in (11. ip given by (11.41) we have 



L(q) = log^det (Gl^Gij 



\?/2 



■^ d c x i(0,l) 



(2.26) 



Now with Pk the probability density function of the nonzero eigenvalues for matrices XX' , 
with X an element of S~ 1 / 2 A/^ xfe (0, 1) as in (I2.17I) - (I2.19I) . this can be rewritten 



L(q) = log ( \\ 



(2.27) 



The average in (I2.27P is given by (I2.2ip with k = 1, z = and \i = q/2, telling us that 



L(q) = log T(l + g/2) 



1 



det 



[1/ 



-9/2] 
i-ll 



3 = l,...,d 



(2.28) 



In the case that E = J d , (12T26|) is a function of G^G^, so a change of variables analogous 
to that used in the derivation of ( \2.7h gives 



L(q) = log^-V^o = log V{q ' 2 ( ^ d) 

T(d) 



(2.29) 



3 Lyapunov exponents for diffusing complex matrices 

The Lyapunov exponents as defined below (11.21) relate to the dynamics of the linear system 
specified by the difference equation Xi+\ = Ai + \Xi for given xq. As emphasized in [19], the 
continuos counterpart of this setting is the matrix stochastic differential equation 

dX(t) = AX(t)dt + BX(t)dW{t) 

where A and B are fixed dxd matrices, and W(t) is a d x d matrix with complex Brownian 
entries. 

In the case that A and B commute, one has 

X(t) = exp ((A - ^B 2 )t + BW(tfjX(0). 

Of particular interest is the case A = \B 2 , B = I so that X(t) = exp(W(t))X(0). Suppose 
furthermore that W(t) = W\(t) + iW2(t) with Wi(t) and Hermitian matrices of 

complex Brownian motions. To specify the latter, let GUEd(0, a) denote the probability 
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density on d x d Hermitian matrices H proportional to exp(— Tr H 2 /a 2 ). We then require 
that W}(1) has probability density proportional to GUE rf (0, aj), for j = 1, 2. An analogous 
specification of W(t) has been given in [TH] for W(t) consisting of real Brownian entries, and 
decomposed as W(t) = Si(t) + 52 (t), where Si(t) is symmetric, and 52 (t) antisymmetric. 
Generally the matrix exp W (1) can be constructed as 

e W(l) = lim e C m (m) e C m (m-l) . . . e C m (l) 
m— >oo 

where C m (j) := W(j/m) — W((j — l)/m). In the above specification of C m (j) is 

independent of j and distributed as C/m 1 / 2 , with C specified as in f ll.lOp . Consequently 
expW(l) has the distribution of Ai for as specified by fll.lOp . 

With this fact established, the argument leading to ( 12. 5 p can be used to show that for 
expW(l) 

Hi + ■ ■ ■ +fi k = lim m(logdet (Ej xk e c ' ^e c ^E dxk )) , 

where the average is over matrices C = H± + iH 2 , with Hj G G UE^ (0, crj) (j = 1,2). 
Straightforward expansion in powers of l/\/m reduces the RHS to 

i(Tr (j5& fc (C* + C) 2 £ rfxfc - (£j xfe (Ct + C)i? dxfc ) 2 ) ) 

= KE £ w) = »fE(<i -■»•)■ 

i=l j=fc+l j'=l 

Consequently we have 

fji k = of(d-2k+l) (3.30) 

(cf. [191 e Q- (15)]). Note that this is independent of 02, and that each Lyapunov exponent 
vanishes for <j\ = 0, corresponding to C — %Hi. This latter point follows from exp W(l) 
then being a diffusion on U(d), and so the modulus of the vectors is unchanged under the 
corresponding flow. 



4 Discussion 

Consider the case k = 1 of (12. ip and thus the maximal Lyapunov exponent. According to 
Proposition [21 for Aj, given by (II. 4p . the exact value of the maximal Lyapunov exponent is 



1 



/'■i 



2 Ui< t<j < d (yj - vd 



det 



[log2/i]i=i,...,<i 

IvT 1 



i=2,...,d 
j=l,...,d 



2 7 



d 

-5E 



2 ^ TT ffl 
j=i LLi=i 



;7, 



(4.1) 
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where 7 denotes Euler's constant. In obtaining the first line the fact that ^(1) = —7 has 
been used, while the second line follows from the first by expanding the determinant by the 
first row and using the Vandermonde determinant formula. Note that replacing 1— > <J~ 2 yi 
for each i — 1, . . . , d changes ji\ by ^\ \-t /ii + logo". To see this from the first line in (14. ip 
requires using the Vandermonde determinant formula, while in the second line one requires 
the identity 

d 

E — — = 1 

j= i Ui=i^ 3 (i-yj/yi) 

(sec e.g. [12, displayed equation below (4.153)]). This change to \i\ is consistent with the 
corresponding mapping of the matrices Ai 1— > aAi. 

For given distinct {yj} we can use the second equation in (14.11) to give a numerical value 
of 71. For example, with d = 2, y\ = 1, 2/2 = 1/4 we obtain [ii = | log 2— 7/2 = 0.63558 
In this case 

1 
2 



Ai 



G1 (4.2) 



where G\ is a 2 x 2 complex Gaussian matrix with entries standard complex normals. 
With xo = [1 0] T , let us define cc, = Aiy^i (i = 1,2,...), where $ := Xi/\xi\. For a given 
realization of {Ai}, it follows from (12. ip with k = 1 that 

fii = lim — \ logic?,!. (4.3) 
i=i 

Moreover, straightforward working shows that — Y2T=i l^il ^ 18bS a Gaussian distribution 



m 1—13 

with mean /ii and a standard deviation a proportional to 1/m 



1/2 



a^-^flogdetfe^) 172 ) 2 ) -ai?) 
m\\\ \ J J I A^ c xl (o,i) / 

This then provides a simple to implement Monte Carlo estimation of \i\ [10]. In the present 
setting, with m = 10 6 we obtained the estimation fi\ w 0.6341. 

Below (II. 8p it was commented that for the average value of the sum of Lyapunov 
exponents to have a well defined limit for d — > 00 it was necessary that the eigenvalues 
{y m } have the scaling form y m /d = Y(y m /d) with Y(0) = 0. Under this circumstance it 
is well known (see e.g. [221 Th. 7.2.2]) that the eigenvalue distribution of G ct EG c tends 
to a well defined nonrandom limit with density ity(t) say, supported on some interval 
I C 1R + . According to a result of Newman [18] (see also [H]) , one then has linid-s.oo e Ml = 
(f T tuyitydt) 1 / 2 . To derive this from (14. ip does not seem possible, although for a given Y(x) 
(14. ip can be used to give a numerical estimation of fi\. For example, with Y(x) = 1 + x, 
computation of (14. ip with d = 5000 (using high precision arithmetic to avoid catastrophic 
cancellations) gives fix ~ —0.183. 

The case k = d of (12. ip . corresponding to the smallest Lyapunov exponent, admits a 
form very similar to the second expression in ( 14. ip . Thus expanding the determinant by 
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the final row and using the Vandermonde determinant evaluation gives 



2 ^ 1X1=1,1^(1 - yj/j/i) 2 

note the interchange of the indices in the denominator relative to (14. ip . With d — 2, yi — 1, 
y 2 = 1/4 as corresponds to (I4.2p . this gives /12 = — § log2 + |(1 — 7) = —0.019656 . . . , and 
thus Hi + ji2 = log 2 + 1/2 — 7 = 0.615932 . . . (note that this is consistent with (jl.8p ). 
Let 4 1} = [1 0] T and xf = [0 1] T , and define x[ P) = A iVf\ (p = 1,2, i = 1,2,...) 
where {$ ,Vi} is obtained from {x^, xf}{\ by the Gram-Schmidt orthonormalization 
procedure. For a given realization of {A}, the analogue of ( 14. 3 j) is then 



/i 1 +/i 2 = lim - Vlogdet^F,) 1 / 2 



where lj is the 2x2 matrix with columns given by and This formula without 



the limit suggests a Monte Carlo estimation of /ii + 1x2 [IQ]. In the present setting, with 
Ai given by (14. 2p . and choosing m = 10 6 gave \i\ + ~ 0.6146. 
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